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Abstract — Timing  synchronization  is  a  preeminent  challenge 
in  ultra-wideband  impulse  radios  (UWB-IRs).  The  conventional 
all-digital  synchronization  methods  encounter  some  formidable 
implementation  difficulties  such  as  high  rate  sampling  and  high 
complexity  RAKE  structure.  To  avoid  these  challenges,  semi¬ 
analog  methods  have  been  motivated  recently.  We  have  recently 
proposed  a  code-assisted  blind  synchronization  (CABS)  algorithm 
to  realize  timing  synchronization  blindly  with  the  help  of  the 
discriminative  property  of  both  time  hopping  codes  and  well- 
designed  polarity  codes.  The  algorithm  requires  sampling  at  the 
frame  rate  only  and  bypasses  channel  estimation  during  the 
synchronization  phase.  This  paper  analyzes  the  identifiability  and 
both  the  probability  of  acquisition  and  the  the  mean  square  error 
performance  of  CABS  analytically.  A  data-aided  code-assisted 
synchronization  (CAS)  algorithm  is  also  proposed  and  a  modified 
version  of  CABS  which  relies  solely  on  the  time  hopping  code  is 
investigated. 

I.  Introduction 

Recently,  ultra- wideband  (UWB)  radio  has  been  attracting 
increasing  interest  due  to  its  unique  advantages  such  as  precise 
ranging  resolution,  enhanced  obstacle  penetration  capability, 
and  low  cost  transceiver  hardware  (see  [1]  and  references 
therein).  All  these  unique  advantages  make  UWB  radio  a 
good  candidate  for  wireless  applications  such  as  location- 
aware  communications  and  wireless  sensor  networks. 

To  realize  the  unique  theoretical  advantages  of  UWB,  many 
physical  and  MAC  layer  problems  have  to  be  overcome. 
Timing  synchronization  is  one  of  the  most  crucial  issues 
since  its  accuracy  has  a  great  impact  on  system  throughput, 
capacity  and  bit  error  rate  [2] [3] [4] [5].  Indeed,  due  to  the 
transmission  of  low  power  and  ultra-short  duration  pulses, 
timing  synchronization  is  very  challenging  in  UWB  impulse 
radios  (IRs),  especially  given  the  low  complexity-low  cost 
philosophy  of  UWB.  Synchronization  algorithms  following  an 
all-digital  route  have  to  cope  with  some  prohibitive  implemen¬ 
tation  hurdles  such  as  the  requirement  of  an  extremely  high 
speed  analog-to-digital  converter,  estimation  of  a  large  number 
of  channel  parameters,  and  a  RAKE  receiver  structure  with  a 
large  number  of  fingers  [6].  This  motivates  the  research  of 
designing  low  complexity  semi-analog  methods[7][8].  In  [8] 
we  proposed  a  code-assisted  blind  synchronization  (CABS) 
algorithm  based  on  the  discriminative  property  of  a  well- 
designed  polarity  code  and  the  time  hopping  (TH)  code  of  the 
desired  user.  With  the  discriminative  nature  of  both  codes,  the 


algorithm  can  achieve  fast  and  accurate  synchronization  even 
in  the  presence  of  multi-user  interference,  while  avoiding  high 
rate  sampling  and  bypassing  challenging  channel  estimation. 

The  algorithm  was  validated  by  both  intuitive  explanation 
and  numerical  simulations  in  [8].  In  this  paper,  we  extend 
our  work  by  providing  an  extensive  analytical  performance 
analysis.  Identifiability  is  first  analyzed.  Then,  for  coarse 
symbol-level  synchronization,  we  study  the  probability  of 
acquisition,  and  a  lower  bound  on  this  probability  is  derived 
in  closed-form.  The  performance  of  fine  synchronization  is 
assessed  using  the  mean  squared  error  (MSE),  and  an  upper 
bound  on  the  asymptotic  MSE  (AMSE)  is  derived.  We  also 
propose  a  data-aided  version  of  CABS,  called  code-assisted 
synchronization  (CAS)  that  provides  better  performance  than 
CABS.  To  relieve  the  challenging  polarity  code  design  in 
CABS,  we  investigate  a  modified  version  which  relies  on  the 
TH  code  only,  and  compare  its  performance  with  the  initial 
CABS  algorithm  by  simulations. 

This  paper  is  organized  as  follows.  Section  II  defines  the 
models  that  describe  the  UWB  signal  and  the  propagation 
channel.  In  Section  III  we  review  the  CABS  algorithm  and 
introduce  the  CAS  algorithm.  The  identifiability  of  CABS  is 
investigated  in  Section  IV,  followed  by  the  analytical  perfor¬ 
mance  analysis  in  Section  V.  Numerical  simulation  results  and 
discussions  are  given  in  Section  VI.  Finally  we  conclude  our 
work  in  Section  VII. 

Notation'.  [_-J  stands  for  integer  floor  operation;  ( -)b  denotes 
the  modulo  operator  with  base  B,  E{-}  and  var{-}  represent 
the  statistical  expectation  and  variance  operators,  respectively. 

II.  System  Model 

In  typical  UWB-IRs,  every  information  symbol  is  conveyed 
by  a  sequence  of  Nf  pulses,  p(t),  each  with  ultra  short 
duration,  Tp,  of  the  order  of  nanosecond.  During  each  frame 
of  duration  Tf,  one  data  modulated  pulse  is  transmitted, 
resulting  in  a  symbol  duration  of  Ts  =  NfTf.  The  two 
most  popular  modulation  schemes  used  in  UWB-IRs  are 
binary  pulse  amplitude  modulation  (BPAM)  and  binary  pulse 
position  modulation  (BPPM),  which  were  compared  in  terms 
of  performance  and  implementation  complexity  in  [9]  [10], 
Since  BPAM  outperforms  BPPM  in  terms  of  bit  error  rate 
(BER),  we  focus  on  BPAM  UWB  here,  and  the  extension 
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of  the  results  to  BPPM  UWB  is  straightforward.  In  order  to 
smoothen  the  spectrum  as  well  as  to  enable  multiple  access, 
periodic  pseudo-random  TH  codes  are  used  to  time-shift  the 
positions  of  the  pulses  by  integer  multiples  of  chip  duration 
Tc  :=  [jf\,  with  Nc  denoting  the  number  of  chips  per 
frame  within  the  associated  frames.  TH  was  shown  to  be  a 
good  spectrum  smoothing  and  multiple  access  scheme  [11], 
especially  for  the  low  data  rate  applications. 

Therefore,  in  the  conventional  TH-UWB  radios  the  pulses 
across  the  frames  of  each  symbol  have  the  same  polarity.  As 
pointed  out  in  [8],  the  polarities  of  the  pulses  within  each 
symbol  can  also  be  modified  by  making  use  of  well-designed 
polarity  codes  to  ease  synchronization.  Thus  the  transmitted 
signal  in  a  point  to  point  link  can  be  expressed  as 

-too  Nf  —  1 

u(t)  =  ^2  ak  22  diP “  kTs  “  iTf  _  CiTc )  (1) 

fc= 0  2—0 


III.  Review  of  CABS  Algorithm 

CABS  algorithm  aims  at  estimating  blindly  the  signal’s 
propagation  delay,  to,  without  channel  state  information.  Blind 
synchronization  algorithms  are  desired  in  many  potential 
UWB  radio  applications  such  as  wireless  sensor  and  ad  hoc 
networks  where  training  may  not  be  available.  The  complexity 
of  estimating  a  large  number  of  parameters  of  a  typical  UWB 
channel  is  prohibitive,  thus  bypassing  channel  estimation  dur¬ 
ing  the  synchronization  phase  is  desirable.  The  idea  behind 
CABS  is  to  rely  on  the  discriminative  property  of  both  the 
polarity  code  and  the  TH  code  to  adjust  the  candidate  time 
shift  until  a  maximum  received  signal  energy  is  captured. 

Since  blind  methods  can  only  provide  synchronization  up 
to  a  symbol  duration,  we  restrict  t0  in  the  interval  [0,  Ts). 
Then,  the  CABS  algorithm  can  be  expressed  as  the  following 
optimization 

fo  =  arg  max  J(t)  (3) 

re[0,Ts) 


where  {a^l’s  are  the  BPAM  information  symbols  taking 
values  ±1  with  equal  probability,  {di}fjo  1  is  the  bipolar 
representation  of  a  binary  code  (i.e.  d.t  =  ±1),  and  {c,  }^jf0  1  is 
the  TH  sequence  whose  elements  are  integer  values  randomly 
chosen  from  [0,-ZV/,  —  1],  where  Nh  <  Nc.  Note  that  if 
d,  =  1,  Mi,  u(t)  reduces  to  the  usual  transmitted  TH  BPAM- 
UWB  IR  signal. 

The  signal  u(t)  propagates  through  an  L-path  fading  chan¬ 
nel  with  {7;}  and  {t/}  representing  the  amplitude  and  the 
delay  of  the  /th  path,  and  without  loss  of  generality  we  assume 
To  <  Ti  <  ...  <  tl-i-  The  channel  is  also  assumed  to  be 
quasi-stationary,  i.e.  {7 j,  T;}[l(J  are  invariant  over  a  burst  of 
several  symbols.  Therefore,  the  channel  impulse  response  can 
be  denoted  by  h(t)  =  7 ~  Ti)-  The  targeted  timing 

information  refers  to  the  delay  of  the  first  arriving  path  tq. 
Thus,  we  isolate  To  from  the  channel  impulse  response  and 
define  the  dispersed  pulse  as  gf)  :=  7 ~  Ti, 0)  with 

T[o  :=  t;  —  To,  Ml  £  [0,  L— 1]  denoting  the  relative  delay  of  the 
/th  path  to  the  first  path.  Therefore,  the  maximum  dispersion 
of  the  pulse  is  Tg  =  Tp  +  tl- i,o- 

The  signal  at  the  output  of  the  receive  antenna  can  then  be 
written  as 


+00  Nf  —  l 

r(t)  =  Y  22  dig "  kTs  ~  iTf  ~  CiTc  ~  T°)  +w(t) 

k—0  i= 0 

s - v - ' 

s(t) 

(2) 

where  s(t)  represents  the  useful  signal  component  and  w(t) 
accounts  for  the  thermal  noise  which  can  be  approximated  as 
a  wide  sense  stationary  white  Gaussian  process  with  double¬ 
sided  power  spectral  density  No/ 2. 

In  this  paper,  we  consider  dense  multipath  propagation 
scenarios  in  which  the  channel  paths  are  closely  spaced.  For 
analytical  tractability  we  also  assume  Tf  >  Tg  +  NhTc  i.e.  no 
inter-frame-interference  (IFI),  although  CABS  was  shown  by 
simulations  in  [8]  to  be  robust  against  IFI. 


with 

J(t) 


1 

K 


dj9r{t)r{t  +  tk'j  +  t) 


2 


dt 


where  tk',j  '■=  k'Ts+jTf  +  CjTc,  K  is  the  number  of  symbol- 
long  segments  of  the  received  waveform  used  for  synchro¬ 
nization,  djgr(t),t  £  [0,  Tf),  with  an  effective  support  smaller 
than  Tf,  is  the  correlation  template  for  the  y  th  frame,  and  Tj  is 
the  integration  interval.  It  is  obvious  that  the  shape  of  djgr(t) 
and  the  value  of  Tj  affect  the  amount  of  the  signal  energy 
and  the  noise  energy  which  are  captured  by  each  integration 
output,  and  thus  affect  the  estimation  accuracy.  For  simplicity 
of  both  implementation  and  analysis,  we  choose  djgr(t)  as 
a  non-return-to-zero  gate  (NRZG)  with  polarity  decided  by 
dj,  i.e.  djgr(t)  =  dj,t  £  [0,7/].  The  optimization  of  T 7  is 
dependent  on  the  actual  channel  and  the  operating  signal-to- 
noise  ratio  (SNR)  value,  and  thus  analytically  untractable.  If 
the  receiver  has  the  knowledge  of  the  exact  maximum  channel 
delay  spread,  Tj  can  be  set  as  Tg.  Otherwise,  Tj  can  be  set  to 
an  upper  bound  of  the  maximum  pulse  dispersion  Tg,  or  even 
to  Tf  -  NhTc. 

It  is  worth  pointing  out  that  if  training  is  allowed,  a 
data-aided  version  of  CABS  (a.k.a.  CAS),  which  can  further 
improve  the  performance,  is  also  feasible.  By  transmitting  a 
training  sequence  with  pattern  of  ak  =  (— l)fc,  CAS  can  be 
formulated  as 


& 

II 

=  arg  max  Jcas{t ) 

re[0,Ts) 

fTI 

Nf-1 

Jcas{t)  =  \ 

djgr(t)f(t  +  jTf  +  CjTc) 

J  0 

3= 0 

where  f(t)  :=  ^  X^fe/=i(-l)fcV(f  +  fcTs  +  t),  t  £  [0,TS] 
is  a  symbol-long  segment  of  the  averaged  received  waveform. 
Although  analytical  analysis  for  CAS  is  possible,  in  this  paper 
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we  only  focus  on  CABS;  extensions  of  these  results  to  CAS 
can  be  obtained  using  the  same  reasoning. 

IV.  IDENTIFI ABILITY  ANALYSIS  FOR  CABS 

In  [8],  the  validity  of  CABS  was  verified  via  an  intuitive 
illustration  and  numerical  simulations,  and  in  this  section 
we  provide  an  analytical  verification  of  the  identifiability  of 
CABS.  Since  the  analysis  when  the  pseudo  random  TH  code 
is  included  is  untractable,  to  simplify  the  noise-free  part  of 
the  objective  function,  we  assume  {cj}f,0  1.  For  analytical 
simplicity,  we  also  assume  I)  =  Tg\  the  ambiguity  introduced 
by  choosing  Tj  fz  Tg  will  be  discussed  later  on. 

Therefore,  the  objective  function  given  in  (3)  can  be  split 
into  three  terms  by  substituting  r(f)  by  (2) 

J(t)  =  Js(t)  +  Jni(r)  +  JU2(t),  (5) 


where 


Js(t) 


Nf-1 

Y.  djs(t  +  k'Ts  +  jTf  +  t) 


3=0 


2 


dt, 

(6) 


K  rp  Nf  —  1 

Jn i(t)  =  jfY  [  £  djlS(t  +  k'Ts+jiTf+T) 

k'  =  lJo  ji  — 0 
Nf-1 

x  £  dj2w(t  +  k'Ts  +  j2Tf  +  r)dt  (7) 

J2=0 


and 


Jn2 (T)  — 


1  K 

-y 

K  ^ 

fc'= 1 ' 


Nf-1 

Y,  djiv(t  +  k'Ts  +  jTf  +  r) 

3=0 


dt. 


(8) 

Since  the  time  delay  To  and  the  candidate  time-shift  r  are 
both  within  the  range  of  [0,  Ts),  the  relative  misalignment 
between  them  is  r  —  To  £  (— Ts,  Ts ).  For  clarity  of  exposition, 
we  consider  two  cases:  t  <  tq  and  r  >  To,  separately. 
Under  each  case,  we  show  that  the  noise-free  term  Js(t)  is 
maximized  if  and  only  if  (iff)  t  =  tq,  which  verifies  the 
identifiability  of  CABS  algorithm. 


A.  T  <  To 

In  this  case,  the  relative  misalignment  between  r0  and  r  can 
be  denoted  as  tq  :=  To  —  t,  with  fo  £  [0,  Ts),  and  can  be  split 
into  a  frame-level  misalignment  quantity  and  a  sub-frame-level 
misalignment  quantity  as:  no  =  [fyj,  £  [0,iV/  —  1],  and 
eo  =  To  —  no Tf,  ?o  £  [0  ,Tf).  After  some  manipulations  and 
algebra,  Js(t)  can  be  simplified  to: 


Nf  A 
Ur)  =  AE 


K 


fc'= l 

32 


R2d[-h0\k'}  /  g2{t-e0)dt  (9) 


TRd[~ Ho  —  1;  k']  J  g2 (t  +  Tf  —  eo)dt 


where  Rd[n\k']  :=  ^  Ef=o  1  d3d{j+n)N/  '  If  the 

two  consecutive  BPAM  symbols  dk>  and  ak,+  ^  j+«j  are  the 

same,  Rd[n;  k']  =  [n]  :=  ^  Ei2o_1  djd{j+n)Nf  ■  Other¬ 

wise,  Rd[n-,k']  =  Rj[n]  :=  ^  E^o”1  dJd(j+n)-N  where 

(n)iv  =  (— (n)jv.  Since  n0  €  [0 ,Nf  —  1]  and  e0  £ 
[0,  Tf),  it  is  obvious  that  Js(t|t  <  To)  is  maximized  iff  no  =  0 
and  ?0  =  0,  i.e.  t  =  To,  and  the  maximum  equals  N2Ng  where 
Eg  :=  ffJ  g'2{t)dt  is  the  energy  of  each  channel-dispersed 
pulse. 


B.  T>  T0 


The  relative  misalignment,  in  this  case,  can  be  denoted  as 
t0  :=  t  —  Tq.  Keeping  the  definitions  for  the  other  parameters 
involved  the  same  as  those  in  the  above  case,  and  following  a 
similar  derivation,  we  obtain: 

\R2d[no\k']  f  g2{t  +  e0)dt  (10) 


Js(t)  = 


K 

ifE 

k'= 1 


Nt 


+Rj[ho  +  k']  /  g2(t  —  Tf  +  eo)dt 


and  with  the  same  arguments  made  in  the  first  case  we  can 
conclude  that  Js(t\t  >  To)  reaches  its  maximum  N2Eg  iff 
fio  =  0  and  e0  =  0,  i.e.  t  =  To- 

Combining  these  two  cases  together,  the  CABS  algorithm 
proposed  in  (3)  can  be  corroborated  identifiable  if  Tg  is  known. 
In  practice,  however,  the  exact  value  of  Tg  is  unavailable  at  the 
receiver.  In  this  case,  the  value  of  Tj  can  be  set  as  either  Tg  or 
a  value  which  is  smaller  than  Tg.  If  Tj  >  Tg,  an  uncertainty 
area  of  [—  (T/  —  Tg)  +  To,  To]  exists,  and  if  Tf  <  Tg,  the 
uncertainty  area  depends  on  the  actual  channel  delay  profile. 


V.  Performance  Analysis 


A.  Probability  of  Acquisition 

In  UWB  IRs,  coarse  timing  acquisition  aims  at  finding 
the  boundaries  between  symbols  in  the  received  waveform. 
In  fact,  the  autocorrelation-based  detection  schemes  such  as 
transmitted  reference  (TR)  [12],  and  differentiate  decoding 
(DIFF)  [13]  are  more  sensitive  to  the  frame-level  timing 
mismatch  [14]  than  the  sub-frame  timing  mismatch.  Therefore, 
the  probability  of  acquiring  the  coarse  timing  synchronization 
is  more  of  interest  than  the  exact  MSE  of  the  delay  estimate 
for  autocorrelation  based  receivers.  Here  we  derive  the  lower 
bounds  on  the  probability  of  acquisition  of  CABS  for  the 
coarse  timing  acquisition. 

During  the  coarse  timing  synchronization,  Nf  candidate 
time  offsets  t  =  nTf  with  integer  n  £  [0,  Nf  —  1]  are 
tested,  which  is  the  same  as  setting  the  search  step  TA  =  Tf. 
If  the  synchronizer  can  lock  the  receiver  onto  a  candidate 
offset  that  is  closest  to  the  actual  delay,  we  then  say  that 
coarse  acquisition  is  achieved.  This  coarse  acquisition  problem 
can  then  be  formulated  as  estimating  no  by  the  following 
optimization 

no  =  max  J(nTf), 

n=0...Nf-l 
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so  that  | To  —  noTf\  <  f .  The  probability  of  acquisition  can 
be  defined  as: 

Pa{n0)  =  Pr{n0  =  n0}  (11) 

=  1  -  Pr  |  U  >  JM) 

or 

/+CO  /  j-  fX  o 

/noW  •••  (12) 

-oo  \J  J  —  oo 

N  fn^no  (a:  i ,  •  •  • ,  a:  tv  f — \)dx-[  •  •  •  dxpj  d  dx  o 

where  (J  stands  for  the  union  of  sets,  fno  (x)  is  the  proba¬ 
bility  density  function  (pdf)  of  the  random  variable  J(nfTf) 
and  fnjtno(xi, . . . ,  dxNf-i)  is  the  joint  pdf  of  the  random 
variables  J(nTf),  n  £  [0,  Nf  —  1],  and  n  Y  uq.  Using 
the  central  limit  theorem,  J(nTf),  n  =  0 . . .  Nf  —  1  can 
be  modeled  as  Gaussian  random  variables.  For  analytical 
tractability,  we  assume  they  are  independent  to  each  other, 
although  the  independence  of  these  random  variables  is  not 
actually  valid  due  to  the  long  channel  delay  spread  and  the 
overlap  among  the  observation  windows.  We  also  assume  that 
eo  =  0  to  simplify  the  derivation. 

The  first  lower  bound  on  Pa  is  derived  based  on  Bonferroni 
Inequality  [15]  and  is  referred  to  as  the  ‘Union  Lower  Bound’ 
(P_ub)-  Denoting  Yn  :=  J(nTf )  —  J(n0Tf),n  f  n0,  and 
with  the  Gaussian  assumption  on  J(nTf),  the  Yn’s  are  also 
Gaussian  random  variables  with  mean  myn  and  variance 
cry  )•  The  Bonferroni  Inequality  states  that  for  a  countable 
set  of  events:  Ai, . . . ,  An,  Pr  {[_}■  A;}  <  JT  Pr  {Ai}  where 
the  equality  holds  iff  these  events  are  mutually  exclusive 
(or  independent).  Substituting  Yn  into  (11)  and  invoking  the 
Bonferroni  Inequality,  we  obtain: 


and  |  Rf  [j]  |  =  \Rf  \j  +  Nf]  |  and  thus  in  (13)  the  summation 
over  Vn  £  [0,_/Vy  —  1]  and  n  ?r0  is  equivalent  to  the 
summation  over  —no  £  [1,  Nf  —  1]. 

Since  it  is  well  known  that  this  PUB  derived  based  on  the 
Bonferroni  Inequality  is  very  loose,  we  resort  to  the  Jensen’s 
Inequality  instead,  and  derive  another  lower  bound  which  is 
named  Jensen  Lower  Bound  (PJB).  With  the  assumption  of 
independence,  the  (Nf  —  1)—  fold  joint  probability  distribution 
in  (12)  can  be  written  as  the  product  of  the  probability 
distribution  of  each  random  variable,  and  (23)  can  be  modified 
as 


Applying  the  Jensen’s  Inequality,  we  obtain  that 

r+oo  f*t-l  i-x  o 

o)  /  /no  (*^o)  1 1  /  fn^n^dXn  dXQ 

J-oo  n=0  J- oo 


> 


\  n=u 
Vn^n  o 

-J  -  O+OO  fX  o 

II  /  fn  o  (X0)  /  fn(xn)dxndx  0 

_ n  J  — OO  J—  OO 


Nf~'  o+o 


n— 0 
n/n0 

n  q 

—  no6[l,JV/  — 1] 


my 


—  PjB- 


(16) 


B.  Asymptotic  MSE 

Since  J(r)  is  a  nonlinear  function  of  r,  the  finite-sample 
MSE  analysis  of  f0  is  extremely  complicated.  However,  with 
sufficiently  large  SNR  and/or  I\,  a  closed  form  upper  bound 
on  the  MSE  can  be  derived.  Due  to  the  page  limitation,  the 
details  of  the  derivation  are  omitted  here  (see  [16]),  and  the 
upper  bound  can  be  expressed  as 


Pa(n0)  >  Pub  ~  1  ^ 

nj^no 


(13) 


where 


Q(x)  =  |  ( 


5  ( 1  —  erf 


(*))• 


with  erf(x)  = 


-A-  J  *  )dt  denoting  the  error  function. 

Since  the  J(nTf)  are  assumed  independent  and  e o  =  0,  the 
mean  and  the  variance  of  Yn  can  be  readily  derived  as: 


myn  =  NfEg 


[+n0]  2  -  R~d  [-no]  2 


-  1 


(14) 


and 

yw" 4  rp2  2 

°Yn  =  -Jff  {-M\2  ^  \Rf  i-n0)\2y  (15) 
(2/v-l)|i?+[-n0]|2  2 

2  aw 

,  |^[-n0]|2  +  2/v  2  ,  2NfTg  4 

+ - ^2 - a v>  +  K 

where  ho  :=  n— no  -  Noticing  that,  due  to  the  periodic  property 
of  R+  [j]  and  Rd  [j],  we  have  that  \R+  [j}\  =  \ R+  [j  +  TV/]  | 


E{f02}  =  Pr(fo  £  [0,  ^])E  {-fg  |fo  £  [0,^]} 

+Pr(f0  £  K,0])E{T02|f0  G  K,0]}  (17) 

<  _ < _ 

-  KN)  -  §g\Tg  -  0} 

where  f0  :=  To  — To,  To  £  [0,  £],  and  £  <Tp  is  a  small  positive 
scalar.  This  reveals  that  the  MSE  of  To  converges  to  zero  as 
K  — »  oo  and/or  SNR  -+  oo.  It  is  worth  pointing  out  that  to 
reach  (17),  or  in  other  words  to  achieve  the  convergence  of 
the  MSE  of  To,  a  specific  shape  of  the  elementary  pulse  p(t) 
is  required.  However,  conventional  pulses  such  as  the  second 
derivative  of  a  Gaussian  pulse  which  is  used  in  the  ensuing 
simulations  also  enable  CABS  to  realize  accurate  estimate  of 
TO- 

VI.  Simulations  and  Discussions 

In  this  section,  we  verify  our  analysis  by  numerical  simu¬ 
lations.  In  all  the  ensuing  simulations,  p(t)  is  chosen  as  the 
second  derivative  of  a  Gaussian  pulse  i.e.  p(t)  =  fEp[l  — 
47t(*~°'5)2]  exp{— 27t(  )2 },  where  tp  =  0.4ns  is  the 

delay  parameter  of  the  pulse.  The  energy  of  this  pulse  is 
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normalized  to  unity,  i.e.  Ep  =  1.  Each  symbol  consists  of 
Nj  =  13  frames  with  frame  duration  Tf  =  60 ns.  Each 
frame  is  composed  of  Nc  =  30  chips  and  the  chip  duration  is 
Tc  =  L]^;J  =  2ns.  The  user  specific  TH  codes  are  generated 
by  randomly  taking  integer  values  from  the  range  of  [0,  N^— 1] 
with  Nh  =  16,  and  the  codes  are  periodic  with  period  Nf.  The 
binary  polarity  code  we  used  is  the  same  as  the  one  we  used 
in  [8].  The  multipath  channel  model  simulated  is  ‘CM1’  from 
the  IEEE  802.15  working  group  [17]  with  maximum  excess 
delay  of  39 ns,  and  normalized  unity  energy.  When  evaluating 
the  normalized  (with  respect  to  Ts2)  MSE  (NMSE),  we  set  the 
interval  of  each  discrete  time  bin  to  be  TA  =  [-^-J  =  4ns 
where  NA  =  15  is  the  number  of  grid  search  per  frame. 
To  evaluate  the  coarse  acquisition  performance  of  CABS, 
the  time  interval  between  two  consecutive  discrete  search  is 
Ta  =  Tf  =  60ns,  i.e.  NA  =  1. 

In  the  above  sections,  we  assumed  the  pulses  are  not 
time-shifted  by  the  TH  code  in  order  to  make  the  analytical 
derivations  tractable.  In  this  section,  we  first  show  by  simu¬ 
lations  that  the  performance  of  CABS  when  TH  is  adopted 
is  consistent  to  that  when  TH  is  absent.  The  NMSE  and  Pa 
of  both  cases  are  plotted  versus  Ep/N0  in  Fig.l  and  Fig. 2, 
respectively.  It  can  be  seen  that  the  two  schemes  have  very 
similar  performance.  Therefore  the  analytical  results  obtained 
can  also  be  used  to  evaluate  the  performance  of  CABS  when 
TH  is  adopted.  In  Fig. 2,  the  lower  bounds  of  Pa  derived  in 
Section  V  are  also  shown.  As  expected,  compared  with  P_UB, 
P_jb  is  tight  especially  at  low  and  medium  SNR  values,  and 
can  be  used  to  analytically  predict  the  acquisition  performance 
of  CABS. 

We  then  compare  the  NMSE  of  CABS  with  that  of  CAS 
in  Fig.  3.  CAS  offers  higher  estimation  accuracy  than  CABS 
especially  at  low  SNR  values,  and  medium  to  large  values 
of  I\ .  Thanks  to  the  special  training  pattern,  extra  averaging 
of  the  analog  signal  is  possible  before  the  integration  for 
CAS,  which  further  eliminates  the  noise.  The  price  paid  for 
achieving  this  performance  gain,  however,  is  some  extra  analog 
delay  lines  and  the  loss  of  data  rate  to  enable  training  sequence 
transmission. 

In  practice,  the  synchronizer  has  no  knowledge  about  the 
starting  time  of  each  frame  during  the  coarse  acquisition  stage, 
and  the  time  positions  at  which  J(r)  is  evaluated  can  be 
anywhere  within  each  frame,  i.e.  ?o  =  0  does  not  hold  true 
in  general.  In  this  case,  the  definition  of  acquiring  the  coarse 
timing  synchronization  given  in  the  previous  sections  needs  to 
be  modified.  In  fact,  if  the  synchronizer  can  lock  the  receiver 
onto  the  time  whose  difference  to  the  actual  start  of  each  sym¬ 
bol  is  smaller  than  one  frame,  i.e.  |fo  —  ro|  <  Tf,  the  coarse 
timing  acquisition  is  considered  successful.  Fig.4  illustrates  the 
probability  of  f0  £  ( To  — 2/ ,  t0+7/)  when  CABS  algorithm  is 
exploited  for  coarse  timing  synchronization  with  eg  £  [0,7/). 

Since  one  of  the  gists  of  CABS  algorithm  is  the  symbol- 
periodic  repetitive  time  position  pattern  of  the  received  wave¬ 
form,  we  are  motivated  to  investigate  a  modified  version  of 
CABS  which  requires  no  need  for  the  polarity  code.  Fig. 5 
and  Fig. 6  show  the  comparison  of  performance  between  the 


E  /N  [dB] 


Fig.  1 :  NMSE  of  CABS  with  TH  vs.  without  TH 


Fig.  2:  Pa  of  CABS  with  TH  vs.  without  TH 

original  and  the  modified  CABS.  Surprisingly,  the  two  have 
almost  identical  NMSE.  The  coarse  acquisition  performance 
of  the  modified  CABS  is  slightly  degraded  by  removing 
the  polarity  code.  This  is  intuitively  explainable.  The  NMSE 
predicts  the  performance  of  the  synchronizer  during  the  fine 
timing  synchronization  stage,  and  the  probability  of  acquisition 
reflects  the  performance  of  the  synchronizer  during  the  coarse 
timing  synchronization  stage.  When  the  coarse  timing  acquisi¬ 
tion  is  achieved,  the  discriminative  effect  of  the  polarity  code 
disappears  during  the  fine  synchronization  phase.  However, 
the  polarity  code  does  affect  the  coarse  timing  performance 
of  CABS.  With  this  exposition,  the  challenge  of  designing 
good  polarity  codes  is  relieved  and  CABS  can  be  enabled 
for  conventional  TH-UWB  systems.  However,  if  the  TH  code 
has  the  property  that  there  exists  any  shifted  version  of  the 
TH  code  which  is  identical  to  the  code  itself,  there  will  be 
ambiguity.  Nevertheless,  the  chance  for  the  TH  code  to  have 
such  property  is  very  small  since  the  value  of  Nf  is  usually 
several  tens  or  hundreds. 

VII.  Conclusions 

In  this  paper,  we  extended  our  work  on  CABS  algorithm 
which  was  initially  proposed  in  [8]  by  providing  an  exten¬ 
sive  analysis.  We  analytically  verified  the  identifiability  of 
CABS,  and  evaluated  the  asymptotic  MSE  and  probability 
of  acquisition  during  the  coarse  timing  synchronization  stage. 
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Fig.  3:  NMSE  comparison  between  CABS  and  CAS 


Fig.  4:  Probability  of  acquisition  with  sub-frame 
misalignment 


A  data-aided  version  of  code-assisted  synchronization  (CAS) 
algorithm  which  outperforms  CABS  was  also  proposed.  To 
relieve  the  polarity  code  design  challenge  and  to  enable  CABS 
in  a  conventional  TH-UWB  system,  we  also  investigated  a 
modified  version  of  CABS  which  requires  no  polarity  code. 
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